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ABSTRACT 

The dissipation of energy from sound waves and weak shocks is one of the most promising 
mechanisms for coupling AGN activity to the surrounding intracluster medium (ICM), and 
so offsetting cooHng in cluster cores. We present a detailed analysis of the weak shock found 
in deep Chandra observations of the Perseus cluster core. A comparison of the spectra either 
side of the shock front shows that they are very similar By performing a deprojection analysis 
of a sector containing the shock, we produce temperature and density profiles across the shock 
front. These show no evidence for a temperature jump coincident with the density jump. To 
understand this result, we model the shock formation using ID hydrodynamic simulations 
including models with thermal conduction and 7 < 5/3 gas. These models do not agree well 
with the data, suggesting that further physics is needed to explain the shock structure. We 
suggest that an interaction between the shock and the Ha filaments could have a significant 
effect on cooling the post-shock gas. 

We also calculate the thermal energy liberated by the weak shock. The total energy in the 
shocked region is about 3.5 times the work needed to inflate the bubbles adiabatically, and the 
power of the shock is around 6 x 10'*^ erg s^^ per bubble, just over 10"*^ erg s^^ in total. 

Key words: X-rays: galaxies — galaxies: clusters: individual: Perseus — intergalactic 
medium — cooling flows 



1 INTRODUCTION 

The baryonic content of galaxy clusters is dominated by a hot 
\ 10^ - 10* K plasma— the IntraCluster Medium (ICM)— which 
■ radiates via bremsstrahlung and line emission processes, making 
' clusters luminous X-ray sources. Observationally the population of 
clusters can naturally be divided into two classes; cool-core clus- 
ters which show a bright central peak in their surface brightness and 
non-cool-core clusters which do not show sharply peaked emission. 
It is believed that the different structure of the two classes is a result 
of different merger histories, with cool-core clusters being more 
relaxed structures which have developed their dense cores without 
significant disruption. 

Measurements of the central luminosities of cool-core clus- 
ters show that they have cooling times much shorter than the age 
of the univ erse. This led to the development of the cooling flow 
model (see lFabiai]ll994l for a review) in which the cooling of gas 
in central regions causes it to lose pressure support creating an in- 
flow of gas toward the centre. The cooled gas was expected to be 
observed both as low temperature X-ray emitting gas and as cold 
gas in the nucleus of the cluster. However it has long been known 
that the star formation rates in central cluster galaxies are below 
the rate expected in the cooling flow picture and, with the launch 
of the latest generation of X-ray satellites, it has been confirmed 
that observed rates of mass deposition from the hot phase are much 
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small er than the expected for a cooling flow jPeterson et alj[200ll 
l2003h . with the core temperature typically falling to about a third 
of the virial temperature. Instead, it is generally accepted that there 
is some heat source counteracting the cooling of the X-ray gas. The 
nature of this heat source is still an open question but the Active 
Galactic Nucleus (AGN) in the cluster centre has emerged as the 
leading contender. 

Observationally, the centres of cool-core clusters often show 
highly disrupted morphologies, in contrast to their overall smooth 
profile. Many cool-core clusters shows "bubbles" — regions with 
depressed X-ray surface brightness profiles indicating that they are 
cavities in the ICM. These cavities are filled with relativistic plasma 
and are evacuated by relativistic jets produced by accretion onto 
the central AGN, thereby coupling the gravitational energy released 
within a few Schwarzschild radii of the central black hole with the 
cool-core region which has a typical diameter of ~ 100 kpc. Since 
the rate of AGN accretion is presumably related to the mass drop 
out rate from the ICM, AGN heating has the potential to be self- 
regulating, thus explaining the ubiquity of cool-core clusters. 

Recent studies of the volumes of observed cavities have found 
that the pAV work done in inflating the cavities is, i n most cases, 
sufficient to offset the X-ray cooling of the cluster gas feirzan et al.l 
2004, 2006; Dunn & Fabian 2006; Rafferty et all l2006h . A natural 
consequence of the expansion of these bubbles is the formation of a 
shock wave which carries the work done in expanding the bubbles 
away from the cluster core. Weak shocks have been observed in 
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severa l nearby clusters, in cluding Perseus jpabian et al1l2006l) and 
Virgo jForman et alj2007h . 

Further evidence of weak shocks is seen in the Perseus clus - 
ter as ripples in the surface brightness jPabian et alj|20033 . l2006h . 
These are interpreted to be low-amplitude sound waves propagat- 
ing through the cluster. That there are more ripples observed than 
bubbles may be an indication that the expansion of the bubble is 
not uniform. 

The dissipation of weak shock w aves at the clu ster centre has 
been i nvestigated by ,Mathews et alj ( [2006 ) and Fuiita & Suzukil 
( l2006h . These authors both use one-dimensional ideal gas models 
to conclude that weak shocks dissipate over a small region com- 
pared to the cooling radius of the cluster and find temperature pro- 
files that are centrally peaked, not centrally decreasing as observed 
in real clusters. Nevertheless, the fact that weak shocks and ripples 
are indeed observed implies that the physics of real cluster cores 
must be more complex than is assumed in these models. 

One of the most surprising c onclusions in the deep study of 
the Perseus Cluster ( iFabian et al .11 20061) was that the weak shock 
in Perseus is isothermal; they found no evidence for the 25 per- 
cent jump in temperature across the shock expected from the ob- 
served density jump. Moreover, multi-temperature fits showed pro- 
portionally more cool (2keV) gas inside the shock than outside, 
strongly suggesting that physics beyond that of a simple ideal gas 
is needed to understand the structure of the shock and hence en- 
ergy dissipation in cluster cores. Here we present a m ore detailed 
study of the data first presented in iFabian et al.l j2006l) , concentrat- 
ing on the region around the inner bubbles in Perseus, including 
modelling of the shock wave evolution. We adopt the same cos- 
mological parameters as in that paper; in particular we assume 
Ho = 71kms"^Mpc-\ 



2 STRUCTURE OF THE SHOCK 

The analysis of the shock front in IFabian et alj j2006l) was based 
on a projected view of the regions inside and outside of the shock. 
Whilst this analysis has the significant advantage that it does not re- 
quire any assumptions about the (unknown) geometry of the cluster, 
the projected gas can have a significant influence on the results and, 
in particular, makes comparison with simulat ion difficult. To ge t 
around this, we have reanalysed the dataset of IFabian et all b006h 
in the sector shown in Fig.[T]performing a deprojection analysis. 

2.1 Surface Brightness Profile 

Fig.|2|shows the projected and deprojected surface brightness pro- 
file in the sector. The sector used to generate the profile here was 
quite wide (80° ). In some places the edge of the shock appears to be 
better defined. We examined several sectors ranging down to 8.9°, 
carefully making sure the sector was aligned along the edge of the 
shock, but the width of the shock appeared to remain consistent. 

2.2 Temperature Profile 

Fig. [3] shows the temperature and density profiles deprojected 
in O.larc min (2.2 kpc) annuli, using the spectral deprojection 
method of ISanderTfe Fabian (2007). This is a simple method to 
calculate deprojected count spectra from observed spectra. Work- 
ing from the outside of the cluster, deprojected count spectra are 
calculated by subtracting from the observed projected spectrum the 




Figure 1. The sector of the Perseus cluster used in the deprojection analysis 
shown on the Chandra surface brightness image of the cluster core. The sec- 
tor is chosen to exclude the regions affected by cooler projected emission. 
The centre of the sector is positioned so that the shock front is at constant 
radius, placing the centre of the sector at the centre of the bubble. 

contribution from those shells which lie between the emitting re- 
gion and the observer. Further details of the deprojection method, 
including a comparison with other deprojection codes such as Pro- 
jCT will be provided in a future paper (Sanders, et al., in prep.). For 
our analysis the spectra have been grouped to a minimum of 100 
counts per spectral bin in the outside annulus, with the others hav- 
ing the same binning. Spectral fits are performed by minimising the 
statistic, with the temperature, column density, abundance and 
normalisation held free, and the redshift frozen at 0.0183. 

At the shock front, the density jumps by a factor 1.31 ± 0.04. 
The Mach number and expected temperature jump can then be de- 
rived from the Rankine-Hugoniot relations, assuming the curvature 
of the shock front has a negligible effect. With 7 = 5/3, appro- 
priate if the ICM is gas (as opposed to cosmic ray) dominated, the 
Mach number of the shock is 1.21 ± 0.03 and the expected temper- 
ature jump a factor 1.27 ± 0.03. 

In the deprojected temperature profile, it is apparent that there 
is no significant temperature jump between the first preshock an- 
nulus and the first postshock annulus. The temperature ratio is 
1.03 ± 0.06. This is clearly inconsistent with the prediction of a 
purely hydrodynamic model. 

Although there is no temperature jump at the shock, there is 
some evidence for a temperature peak ~ 3 kpc behind the shock 
front. The ratio of this maximum in the postshock temperature and 
the minimum of the preshock temperature is 1.14 ± 0.07, which is 
marginally consistent with either a flat profile or a standard shock 
jump. However, interpreting this feature as a standard shock re- 
quires some process beyond simple hydrodynamics to explain the 
offset between the shock face and the temperature jump. 

The annulus centred around 11.8 kpc, shows a drop in both 
the temperature and the density, corresponding to a drop in the 
gas pressure. The origin of such a drop is not well understood; 
it suggests an additional non-thermal contribution to the pressure 
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Figure 2. Projected (upper curve) and deprojected (lower curve) surface 
brightness profile in the sector of the Perseus cluster in Fig.[T] The units are 
arbitrary. 
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Figure 3. Electron density (top) and projected and deprojected (circles) 
temperature profiles (bottom) in the NE sector of Perseus. Radii are mea- 
sured relative to the centre of the sector The dotted line shows the position 
of the shock. 
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Figure 4. Comparison of the projected (diamonds)and deprojected (circles) 
spectra in the immediately inside (filled red points)and immediately outside 
the shock (unfilled blue points). The projected spectra are presented with- 
out any scaling, whilst the deprojected counts (circles) are scaled to units 



of 10" 



' keV, to take into account the area of the 



sector. The ratio plot (bottom) gives the ratio of the counts per unit energy 
inside the shock to those outside, scaled to have the same overall normali- 
sation. The projected spectra inside and outside the shock show remarkable 
consistency. 



from e.g. magnetic fields, cosmic rays, or ram pressure in this re- 
gion. Taking only points outside this, the overall temperature profile 
shows evidence for some structure with a constant model providing 
a poor fit (reduced = 2.7). 

Accounting for the need to fit both the actual structure of the 
shock front, and the fact that the temperature profile far in front and 
behind the shock is not well modelled, we choose to proceed by 
comparing models profiles to the temperature points in the region, 
16 — 25 kpc. This encompasses two annuli inside and three annuli 
outside the shock. 



2.2.1 Checks on the Deprojection Analysis 

In order to check that the low apparent temperature jump behind 
the shock is not an artifact of the deprojection procedure, we have 
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Figure 5. Fitted temperature for NEI models with rie 
4.5 keV at various times after tlie sliock has passed. 
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compared the spectra from the annuli either side of the shock front. 
Fig. |4] shows the projected and deprojected spectra from regions 
immediately inside and outside the shock front. The spectra on each 
side of the shock are remarkably similar, both in projection and 
in deprojection, consistent with the small temperature jump found 
from the full deprojection analysis. 

An additional test we made was to use the spectrum from out- 
side the shock as a background for the shocked region. This is 
almost certainly an overestimate for the amount of emission pro- 
jected onto the shocked region. Using this background, the mea- 
sured temperature was still compatible with our deprojected results. 
It appears not to matter how much projected emission is subtracted 
from the shocked region as the spectra are almost identical. 



3 MODELLING AGN-INDUCED WEAK SHOCKS 
3.1 Equations of fluid dynamics 

For the purposes of our model, we assumed the dynamics of the 
ICM are those of a single phase, inviscid fluid in which magnetic 
fields are too weak to be dynamically important and where the 
electron-ion coupling time is much shorter than any timescale of in- 
terest. In this case the equations descri bing the evolution of the IC M 
are those of ideal hydrodynamics (e.g. lLandau & Lifshitjl959h : 



Du ^ 



De 



D 



Dt 



= — V ■ 
P 



(1) 



(2) 



(3) 



where D = ^ -I- it • V is the total (Lagrangian) derivative, u is 
the fluid velocity, p the pressure, p the density, e the specific inter- 
nal energy, T the temperature, and k the conductivity. Fo r a fully 
ionised hydrogen plasma, the conductivity k is given by ( jSpitzeil 
[l95l): 



: 640e5T 



ZlnA 



(4) 



where Z is the charge of the ion, A is the Coulomb logarithm with 
In A ~ 37 in the conditions of interest, e ~ 0.4 accounts for the 
reduction of conductivity in a plasma due to the electric field and 
5t ~ 0.225 accounts for the finite proton mass. The remaining 
symbols have their usual physical meanings. 

In the case of a steep temperature gradient and strong con- 
ductivity, the rate of heat flow predicted from the bulk description 
above can exceed the maximum rate at which the electrons are able 
to transport energy. In this case it is appropriate to replace the con- 
duction term in the energy equation Q by o ne appropriate for sat- 
urated conductivity jCowie & McKeelll977h : 
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2.2.2 Non-Equilibrium lonisation 

The models fitted above assume that the ions in each zone are 
in their equlibrium ionisation state. Following the rapid tempera- 
ture change induced by the passage of a shock, there is a finite 
time before ionisation equlibrium is re-established. During this 
time, the temperature inferred from X-ray spectral fitting will be 
lower than the actual gas temperature, suggesting a possible expla- 
nation for the offset of the temperature ju mp at the shock front . 
This was tested using t he XSPEC NEI model ([Hamilton et al.|T98^ : 
iBorkowski erZlll994 iLiedahl et al.l ll995l ; ISorkowski etail2001l) 
with a PHABS absorbed component to simulate Chandra observa- 
tions at various times after the shock passage and fitting the sim- 
ulated data with the EQUIL model. The temperature fitted to these 
models is shown in Fig. |5] clearly ait = l(f yr the fitted temper- 
ature is equal to the input temperature. Since the shock velocity is 
approximately 1 kpc/10*' yr non-equlibrium ionisation effects can 
account for an offset of no-more than 1 kpc, indicating that non- 
ionisation equlibrium alone is not enough to explain the apparent 
density/temperature offset at the shock front. 



3.2 Shock Propagation in the ICM 

In general, propagation of a shock wave in a highly ionised plasma 
such as the ICM differs from that in an ordinary fluid because of the 
very different sound speeds for the electrons and ions. At the shock 
front, the ions undergo a temperature and density jump given by the 
normal Rankine-Hugoniot relations for a shock of Mach number 
M: 



7 + 1 



P2, shock _ 
pi, shock (7 - 1) + ll^ 



Tl, 



(27M^ ^ (7 - 1)) (2 + (7 - 1)M^) 
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Charge neutrality requires that the electrons are compressed 
adiabatically at the shock front in the same ratio as the ions, leading 
to a temperature difference: 

7-1 
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P2, shock 
Pi, shock 



(8) 
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The two components are brought into thermal equihbrium be- 
hind the shoc k over a timescale comparab le to their thermal equili- 
bration time dZel'Dovich & RaizeJl967h : 



1.1 X lO^^MA 
(ni/ cm-3)Z2 In A 



keV J 



3/2 



(9) 



Where Tci is the electron-ion equilibration timescale, A the 
atomic mass of the ions, and rii the number density of ions. 
For the case of the Perseus cluster core with rii ~ 0.03 cm~^, 
T ~ 3.5 keV, A = 1.29, ^ = 1.2 and M = 1.2, 
electron-ion equilibration occurs approximately 3.2 kpc behind the 
shock, consistent with the apparent offset between the tempera- 
ture and density jumps at the shock front. In the case of a weak 
shock, however, the heating is largely due to adiabatic compres- 
sion which acts equally on the electron and ion components and so 
from ^ the temperature difference is only ~ 1 percent. Therefore 
electron-ion non-equilibration is unlikely to explain the observed 
density/temperature offset and we can approximate the behaviour 
of the ICM in the shocked region as that of a single-component 
fluid. 



3.3 Conductivity of the ICM 

In general the conductivity of the ICM may differ from the Spitzer 
value of equation Q because even dynamically unimportant mag- 
netic fields are effective in suppressing transport across field lines. 
Therefore the level of conductivity is strongly affected by the ge- 
ometry of the field. In general it is assumed that the field is suf- 
ficiently tangled that, on scales of interest, isotropic conduction 
occurs due to motion of electrons along diverging field lines and 
so the c onductivity is the Sp itzer form reduced by some factor 
/ < 1; iNaravan & MedvedevI l llOOlb derive a suppression factor 
of / = 0.2 based on the assumption of a tangled field. If the field is 
more ordered on smaller scales, the conduction will not be isotropic 
and parallel to the field lines the conductivity will be close to the 
full Spitzer value. 

In the case of the ICM, there is some evidence of ordered field 
structures on scales > a few kpc. Surrounding many central cluster 
galaxies are bright, line emitting neb ulae. In several cases including 
the Perseus (Conselice e t al.ll200lh and Centaurus clusters, long, 
largely radial, filamentary structures (up to 60 kpc in projection), 
are observed which are apparently long lived despite bein g embed- 
ded in the hot ICM (' pabian et alj2003ttlHatch et al.l2005h . In order 
to be stabilised against evaporation by the surrounding gas, thermal 
conduction between the filaments and their surroundings must be 
supp ressed by about an orde r of magnitude more than in a tangled 
field jNipoti & Binnevl2004l) . This is suggestive of a field geometry 
in which the filaments are enclosed in a sheath of radially-coherent 
field lines. 

Increasing radial coherence of the magnetic field toward the 
centre of the cluster is expected if there is radial flow of gas with a 
radially decre asing velocity amplitud e, for example in a (reduced) 
cooling flow ( Soker & Sarazin|[l99Clh . Recent numerical studies of 
bubbles evolution in a magnetised ICM has also shown the mag- 
netic field becomes more ordered in the wake of the rising bubble 
( iRuszkowski et alj|2007bl) . With these results in mind, we assume 
that conduction may proceed radially at up to the full Spitzer rate. 



3.4 Numerical Method 

We have written the fluid equations (|Tll - ([3} in a form appropriate 
to a one dimensional, spherically symmetric system: 
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where q is the appropriate conductive heat flux. These are supple- 
mented with the ideal gas equation of state: 



-kuT 



(13) 



where we have taken the average particle mass /i as a constant set 
to 0.61 appropriate to the chemical composition of the ICM. 

To numerically solve these equations, we used an explicit fi- 
ni te difference method on a La grangian grid, based on the method 
of Ritchmver & Morton] ( 1 1967b and following the implementation 
of Morris ( 2003). An artificial viscosity term was added to the equa- 
tion of motion to correctly model the shock jump conditions. The 
strength of the artificial viscosity was chosen so that the shock was 
smeared over approximately four grid zones. 

For simplicity we neglect the contribution of the cluster galax- 
ies to the gravitational acceleration and assume the contribution 
from gas is a fixed fraction of the dark matter contribution. Assum - 
ing the dark matter follows a NFW profile ( lNavarroetai]|l997h . 
the acceleration may be written: 



ffNFW 



GMviriai log(l + cr/rvinai) - CT / {r^.nji + cr) 
r2 log(l + c) -c/(l + c) 



(14) 



where Afviiiai and r-yiriai are the virial mass and radius respectively 
and c is the cluster concentration. 

The shock wave was generated by moving the position of the 
inner boundary. If the jet supplies energy at a constant rate E then 
it can be shown that the bubble radius varies with time as: 



Rb — 



P 



(15) 



Numerically this formulation is problematic as it requires a 
high (formally infinite) velocity at t = 0. Instead we have adopted a 
simpler model with an initial exponential increase in the expansion 
velocity to a maximimi and a subsequent fall off like cos(t): 



^piston 




< t < tsmooth 
^smooth ^ t <^ tmax 

(16) 



The smoothing time tsmooth has been fixed at 1 Myr; with this 
time, no spurious central heating due to sharp changes in velocity 
are observed, tp is a characteristic piston timescale and ito defines 
the maximum piston velocity. Since we are concerned with only 
a single injection event, we choose tmax to occur when the piston 
velocity first reaches zero, i.e. tmax ~ ^smooth + 



© 0000 RAS, MNRAS 000, 000-000 



6 J. Graham, A.C. Fabian & J.S. Sanders 




50 75 
Radius / kpc 



Figure 6. Deprojected temperature data from the 200ks observation 
iSanders et al. 2004) for a sector of the Perseus cluster containing the shock 
(points) and adopted initial temperature profile (line). The sector used con- 
tains the four regions covering the shock in Fig. 18 of (Sanders et al. 200^. 
This region is similar to that in Fig.[T]but has a significantly greater radial 
extent. The vertical line shows the approximate position of the shock. 
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4 COMPARISON WITH DATA 
4.1 Initial Conditions 

As we are primarily concerned with modelling the observed shock 
in the Perseus cluster, the majority of our models are set up with 
properties close to those of the NE sector of the Perseus cluster 
used in the deprojection analysis. Our initial temperature profile 
was based on the observed temperature in the pre-shock region (i.e. 
r > 20 kpc). To determine the shape of the profile at karge radii , 
we used a fit to the 200 ks data presented in lSanders et alj j2004l) 
covering the shock region, whilst the temperature in the immediate 
vicinity of the shock was ma tched to the data in Fig. [3] The result- 
ing profile, together with the I Sanders et alj ( |2004|) data is shown in 
Fig.m The functional form is: 



Figure 7. Density and temperature after the shock has reached 24 kpc with 
7 = ^ and no conduction (solid line), a Spitzer fraction of 1.0 (dashed 
line) and 2.0 (dotted lines) and deprojected observations in the NE sector 
(points). 

shock approximately reproduced the expected density jump at the 
current shock position. 

When adding conduction to the models, Spitzer fractions, /, 
in the range 0-2 have been considered. Cases with the conductivity 
in excess of the Spitzer value are intended as a simple approxima- 
tion for a situation where the electrons and ions are decoupled. In 
this case only ~ 1/2 of the energy needs to be transported to de- 
crease the observed (electron) temperature compared with the fully 
coupled case. 




The gravitational acceler ation was mode l led us ing the NFW 
component of the profile in iMathews et ( l2006h . which pro- 
vides an acceptable fit to the observed acceleration in the region 
r > 10 kpc which is of interest here. The initial density profile is 
calculated assuming hydrostatic equilibrium and a central density 
rie = 0.05 cm^"^. This gives a gas mass fraction inside 25 kpc of 
~ 6%, so a static gravitational field is expected to be a good ap- 
proximation. The piston timescale tp was taken as 1 x lO'^ yr and 
the smoothing time as 1 x 10® yr. 

The inner boundary was started at r = 5 kpc — approxi- 
mately the distance from the bubble centre to the cluster centre — 
and the piston amplitude adjusted to rpiston = 3 kpc so that the 



4.2 Simulation Results 

For the purposes of comparison with the simulations, we have used 
radii measured from the cluster centre rather than from the centre 
of the deprojected sector, since this more closely matches the setup 
of the simulated cluster. The radius of the shock is then 24 kpc. 

Fig. [7] shows the density and temperature structure of a sim- 
ulation with 7 = 5/3 after the shock has reached the observed 
radius (taken to be 24 kpc). The time taken for this to occur was 
dependent on the level of conduction; in the non conductive case 
the time was ~ 12 Myr, stronger conduction reduces the pressure 
difference across the shock so slowing its propagation. The pis- 
ton amplitude was adjusted so that the observed density jump was 
matched in each case. The addition of conduction with f — 1 and 
/ = 2 reduces the height of the jump to about 65 percent and 50 
percent respectively of the non-conductive case. However conduc- 
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Table 1. Results of numerical simulations for various values of 7, the 
Spitzer fraction / and numbers of injection events. The piston amplitude 
Tp is adjusted so that the density jump at 24 kpc matches the observations. 
The quoted reduced chi-squared is based on a fit of the temperature data 
alone to the five data points in the region 18 — 28 kpc enclosing the shock 
front. 



tion does not appear to change the position of the temperature jump 
relative to the density jump. 

A more quantitative comparison of the data and the models 
can be made by performing a fit of the model temperature profile 
to the data in the region surrounding the shock. To allow for the 
fact that the data shows more structure ahead of the shock than the 
initial conditions of our model, we have selected only those data 
points in the region 18 kpc ^ r ^ 28 kpc and spectrally averaged 
the model temperatures to match the data bins. A simple fit of the 
models to the data products a rather poor fit in each case with the 
reduced ranging from 4. 1 in the case / = 2 to 7.2 in the case 
f — 0. The quality of the fit also depend somewhat on the initial 
temperature profile assumed. For example, maintaining the shape 
of the profile but increasing the central value makes the fits poorer 
in all cases. 

We have also considered models in which the ratio of spe- 
cific heats is less than the standard ideal gas value of 5/3. Lower 
values of 7 may be appropriate if the ICM in the neighbour- 
hood of the shock is significantly contaminated by relativistic 
particles diffusing from the bubbles l lMathews & Brighentill2007l : 
iRuszkowski et al.ll2007ah . In this case the effective adiabatic index 
is given by: 



7eff : 



7Th + ^cr7cr 
1+XcR 



(18) 



with XcR = pcR/pTh jPfronmier et alJl2007h . We consider 7 = 
4/3 as a limiting case of cosmic-ray-pressure dominated gas and 
7 = 3/2, i.e. XcR — 1, as a realistic value for the core of the 
Pe rseus cluster based on th e analysis of the non-thermal emission 
bv lSanders & FabianI JioOTh . 

The temperature profiles for the 7 = 4/3 case are shown in 
Fig. [8] As expected from the standard shock jump conditions, the 
magnitude of the shock in the is substantially decreased compared 
to the 7 = 5/3 case. However conduction appears to show lit- 
tle further effect in this case. The reduced 7 appears to somewhat 
improve the agreement between simulation and data but does not 
produce a statistically acceptable fit in either case. 
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Figure 8. Shock model at r = kpc for the Perseus-like model cluster with 
7=1 and no conduction (solid line), a Spitzer fraction of 1.0 (dashed line) 
and 2.0 (dash-dotted line). 



4.2.1 Models with Multiple Consecutive Outbursts 

By construction our models match the density and the density jump 
at the shock front. Nevertheless, the models considered so far all 
show high density regions inside the shock considerably narrower 
than that observed. 

One way of addressing this issue is to invoke multiple out- 
bursts during a single cycle of AGN activity. Fig.[9]shows the den- 
sity and temperature profiles for set of models in which the single 
outburst has been replaced by two identical outbursts beginning at 
t = Oandt = 4.5 Myr . 

This model successfully reproduces a thicker rim of high den- 
sity gas that extends in approximately as far as the obscuring cool 
gas seen in the data, although it tends to produce a somewhat worse 
fit to the temperature profile than the single-outburst models. 



4.3 The Temperature Offset 

Since our simulations do not convincingly reproduce a tempera- 
ture/density offset at the shock front, we must consider other mech- 
anisms for cooling the gas in the immediate aftermath of the shock. 



4.3.1 Non-Uniform Initial Temperature 

We have assumed a temperature profile that is flat in the core of 
the cluster. However the deprojected temperature profile shows ev- 
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Figure 9. Temperature profile for a model with two episodes of injection 
starting at t = and 4.5 X 10^ with rp = 3 kpc and 2 kpc respectively 
and tp = 1 X 10^ yr in each case. 



idence for oscillations in the temperature of up to ~ 0.5 keV. A 
blob of cold gas producing a temperature fluctuation of around this 
magnitude at the position of the observed density jump can produce 
density and temperature profiles similar to the observed ones — a 
simulation of such a profile is seen in Fig. [TO] However this re- 
lies on fine-tuning of the initial temperature profile, suggesting that 
weak shocks in general would not show such an offset. However, 
the very high resolution observations required to test this are only 
likely to be possible for a very small number of local clusters in the 
foreseeable future. 



4.3.2 Mixing of Filaments 

Thus far, we have ignored the multiphase nature of the Perseus core 
in our analysis. Of particular significance is the reservoir of cool 
gas contained in the large filamentary nebula associated with the 
central cluster galaxy. However, there is observational evidence for 
an interaction between the phases; Fig.ll llshows the X-ray pressure 
and Ha map around the region we have analysed. The filaments in 
this region appear to terminate at the shock front, suggesting the 
cool gas from the filaments is being mixed in with the hot ICM 
in this region. This mixing could have a significant impact on the 
temperature structure of the shock. 

The filaments are known to contain H2 with a tempera- 
ture of 300 — 10000 K and a mass per u nit length of around 
10" - 10^ Mq / kpc jjohnstone et"ai]|2007b. Associated with the 
filaments is soft X-ray emission jpabian et al ]|2003d) , tracing a to- 
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Figure 10. Density (top) and temperature (bottom) profiles at t = (dashed 
lines) and i = 13 Myr (solid lines) for a model with a blob of colder gas 
close to the observed shock radius. 



tal mass of 10 M© within the inner 1.5arcmin (~ 33 kpc) 
iFabian et al.l ( |2006|). There is also a detection of CO in the region 
of the filaments (^Sal ome et al.ll2006h . implying a population of gas 
at 10-lOOK, although it is not clear if this is confined to the fila- 
ments or more dispersed. If it is indeed confined to the filaments it 
implies a mass per unit length of around 10* Mq / kpc. The den- 
sity of the X-ray emitting gas in the region of the shock is approx- 
imately 7 X 10^ Mq / kpc"^ so, assuming complete mixing of the 
behind the shock, a single filament could cool around 1000 kpc^ of 
shock front by 0.5 keV, an area comparable to the surface area of a 
spherical shock centred on the Perseus bubbles and at the observed 
radius. However, in practice the efficiency of mixing is likely to 
be considerably less than 100 percent, so this estimate is an upper 
bound. 

If the post-shock gas is cooled by mixing of filaments, then we 
might expect a spatial variation in the post-shock temperature ac- 
cording to the distribution of filaments along the shock front. Fig. 
[T2I shows the deprojected density and temperature profiles when 
the sector in Fig.[T]is divided into two equal-angle sectors. Outside 
the shock front the two profiles are consistent, but inside the shock 
there is a significant difference in the temperatures for the two sec- 
tors; one is consistent with a flat temperature profile whilst the other 
has a jump of 1.36 ± 0.13 between the preshock minima and post- 
shock maxima. Interestingly, it is the sector closer to the currently 
observed filaments that has the larger temperature jump; this may 
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be indicative of the shock having destroyed previous filaments in 
the cooler sector. 

Assuming the filaments have had a significant effect on the 
temperature structure of the gas, it is not clear why the postshock 
temperature should be roughly the same as the preshock tempera- 
ture; by mixing more or less cold gas into the ICM this could be 
varied from the postshock temperature expected from the Rankine- 
Hugoniot relations to a temperature significantly lower than the 
preshock temperature. A physical mechanism to explain the near- 
isothermality observed would make this model significantly more 
appealing. 



4.4 Repeated Shock Heating 

We now briefly address the issue of r epeated shock heat- 
ing of the cluster co re. Previous authors ^Mathews et al.l [Tood ; 
iFuiita & Suzukill2006h . have found that repeated shocking of the 
cluster core leads to a centrally increasing temperature profile, 
clearly inconsistnet with the observations. This is a result that we 
can reproduce in our model; Fig. [T3] shows the result of repeat- 
edly shocking a cluster, with the initial conditions identical to those 
employed for the models of the Perseus shock, over a period of 
4 Gyr with wave parameters Vp = 2.75 kpc, tp — 2 x 10^ yr and 
fmin = 10 kpc. Clearly the temperature profile produced is not the 
smoothly decreasing temperature profile we expect from observa- 
tions. 

However, we recommend caution in a literal interpretation of 
these results derived from these ID models. There are several lim- 
itations of the models that may cause the discrepancy between the 
model and the data: 

• The ID nature of the models misses several essential features 
of real clusters. The geometry of bubble heating in 3D allows the 
redistribution of material from outside the bubble into the core of 
the cluster and conversely the entrainment of the material behind 
the ri sing bubble draws material away from the centre (e.g. GardlnJ 
|2007|) . Changes in the jet direction and bulk motion of the ICM in 
the cluster potential are also significant effects in ensuring that the 
same cluster gas is not rep eatedly subject to the strongest heating , 
as happens in the ID case (Brtig gen et alj2005l ; lHeinz et al. 200d). 

• The assumption of constant piston amplitude ("bubble size") 
over time is likely to be unrealistic. One would expect a real 
source to release bubbles with a variety of sizes over a variety of 
timescales. Evidence of this in the Perseus cluste r is seen in the spa- 
tial va riation of the ripple amplitude reported bv lSanders & Fabi^ 
( l2007h . In addition to changes in the intrinsic power of the central 
source, the bubble size and release frequency, at constant power, 
will vary with the surrounding gas properties. However the weak- 
ness of this dependence means this is likely to be an insignificant 
effect. 

• Models to date have taken no account of the possibility of the 
mixing of cool gas from the brightest cluster galaxy nebula into 
the ICM. We have shown that this process may plausibly have an 
effect on the temperature structure of the weak shock in Perseus, 
in which case it will alter the temperature structure produced by 
repeated heating episodes. 

• A non-homogeneous population of cosmic rays will affect 
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Figure 13. Density and temperature profiles for a simulation in which a 
Perseus-like cluster is subject to monochromatic wave heating for 4 Gyr. 
The central temperature peak is characteristic of this kind of simulation but 
is inconsistent with the data. 



how the energy dissipated by weak shock waves is distributed spa- 
tially. 

In order to resolve these issues, it is clear that simulations 
with repeated bubble formation events run o ver a period of sev- 
eral G yr in multiple dimensions are required. lOalla Vecchia et aTl 
( |2004|) have run a 3D simulation for 1.5 Gyr in which they are able 
to balance heating and cooling, although their results are dependent 
on the radius over which their bubbles are distributed, with small in- 
jection radii (consistent with observations) leading to a potentially 
problematic increase in energy of the cluster core. However their 
simulation contains neither bulk gas motion not the inflation of the 
bubbles themselves. Sijacki et al. (2007) have also run a cosmolog- 
ical simulation including mechanical feedback from AGN, with the 
energy and radius of the bubbles coupled to the black hole accretion 
rate. They find the injection is able to prevent catastrophic cooling 
in the cluster core, whilst retaining a flat or decreasing temperature 
profile in the centre. 



5 ENERGETICS OF THE SHOCKED REGION 

The amount of mechanical energy imparted to the cool core region 
by the AGN is critical in determining whether the AGN alone can 
offset the cooling flow to the observed levels in cluster cores. If 
the bubble is slowly inflated and thereafter expands adiabatically 
as it rises buoyantly through the cluster core, it releases an en- 
ergy of approximately (7bubbic/(7bubbie - i))pV, i.e. its enthalpy 



l lChurazovetalJl2002h . Therefore studies of AGN heating have 
typically assumed E = 4pK)bservcd, appropriate if the gas inside 
the bubble is purely relativistic. With these assumptions, authors 
have typically found that the available energy is enough to offset 
the cooling in some but not all clusters (e.g. iBirzan et alj j 2004l : 
iDunn & Fabiarill2006l : iRaffertv et alj2006l : iBirzan et alj2006l) . The 
fraction of clusters for which heating is observed to balance cooling 
at the present time is often used to infer a duty cycle for mechanical 
energy input by the AGN. 

Some simulations have suggested that the overall energy re- 
lease is substantially greater than assumed in these s tudies due to 
the contribution from non-adiabatic processes (e.g. Isinnev et al.l 
2007). The high quality data from the Perseus allow us to see the 
high pressure rims around the inner bubbles, allowing an observa- 
tional estimate of the extra energy deposited by irreversible pro- 
cesses during th e bubble expansion . 

Following iDunn et lii ( |2005|) we assume that the Perseus 
bubbles are approximately spherical, of radius rin„cr and that the 
shocked material is confined to a spherical shell of radius rou(cr out- 
side each bubble. The energy dissipated by the shock is then: 



-£'shock 



7g^ 



' Punshocked)'^V" 



(19) 



The value of Pu„shockcd is, of course, an unknown quantity, 
which we must extrapolate from the properties of the cluster out- 
side the shock. In order to minimise the uncertainty associated with 
this extrapolation, we assume that both bubbles are identical and 
that the pressure profile through the whole region surrounding the 
bubble, both inside and outside the shock, is well approximated by 
the profile in that region in the sector we have studied in this pa- 
per. With these assumptions we only need extrapolate the pressure 
profile over the width of the high pressure region. To perform this 
extrapolation, we use a power-law fit to the pressure profile in the 
range 20 — 45 kpc, measured from the bubble centre: 



Punshocked =0.40 



1 kpc 



keVcm 



(20) 



For the shocked pressure we used the values measured in 
the high pressure region r = 10.1 — 19.3 kpc directly (where 
r is measured from the centre of the bubble as in Fig. O. With 
these parameters, the excess energy in the high pressure shell is 
1.6 X 10^^ erg per bubble. If we estimate the power provided by 
the shock as i^shock/tshock where tshock = (rome, - dinner )/Mc, we 
find /'shock ^ 6 X lO^'^ erg s ^ pe r bubble. 

From lDunn & Fabiaij ( |2004|) . the pV work done in expanding 
the bubbles adiabatically are 3.9 x lO'^* erg and 5.3 x 10^** erg 
for the northern and southern bubbles, respectively. Therefore, av- 
eraging over the two bubbles, we find that the energy in the shocked 
region is roughly 3.5 times the pV work done in expanding the bub- 
bles adiabatically, already comparable to the enthaphy of the bub- 
bles. Assuming this result generalises to other clusters, estimates of 
the total AGN-supplied mechanical power available to heat cluster 
cores which assume ApV energy per bubble over their lifetime are 
likely to be reliable. 

The calculated power is somewhat higher than the range 
6 X 10''^ - 2 X 10"* ergs"^ per bubble found bv lPunn & FabiatJ 
(2004) on the basis of the work done expanding the bubbles (i.e. 
taking the energy to be pV, not 4pV) and various estimates of the 
timescale for energy dissipation. It is also comparable to, although. 
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again, slightly larger tha n, the power in the ripples calculated by 
ISanders & Fabianll2007h . 

These estimates have all been ma de on the assumption tha t 
only the thermal pressure is significant. I Sanders & FabianI ( |2007|) 
found that inside 40 kpc Perseus has a significant non-thermal pres- 
sure component contributing about 50 percent of the thermal pres- 
sure. Therefore the energy inside the shocked region may be an 
additional factor '-^ 1.5 higher than derived here. 

It is also possible that excess energy has been channelled 
into other modes such as gravitational energy, kinetic energy or 
turbulence. We expect the contribution of gravitational energy to 
be small; assuming the expansion is approximately symmetric, 
the almost linear nature of the cluster potential over the region 
r < router means that the net change in the gravitational poten- 
tial cancels. Relaxing the condition of symmetrical redistribution, 
the condition of hydrostatic equilibrium gives AP ~ PpAr so 
£^grav/£^thcrniai ~ AP/P ~ 0.1 therefore the gravitational contri- 
bution to the total energy is likely to be small. 



6 DISCUSSION AND CONCLUSIONS 

We have performed a detailed analysis of the weak shock feature 
seen in the core of the Perseus cluster. Deprojection of a sector 
containing the shock shows no temperature rise coincident with 
the density jump. There is a small rise 3 kpc downstream from the 
shock front. Models of shock production in spherically-symmetric, 
ideal gas, cluster atmospheres produce poor fits to the observed 
profile, principally because the observed postshock temperature is 
lowest at the shock front, just where it is expected to be highest. 

The presence of thermal conduction, or a reduced value of the 
adiabatic index from the presence of a population of cosmic rays, 
can produce a temperature jump with very similar magnitude to 
that seen in the data and substantially improve the fit of the shock 
profile to the data. However we did not find a set of parameters for 
which the value of the fit is statistically acceptable. The diffi- 
culty of finding a simple solution to this problem is likely related to 
the fact that the densities and length scales relevant to the Perseus 
shock represent the border between the collisional and collisionless 
regimes. Therefore the non-hydrodynamic physics, such as the in- 
teraction of the electron population with the local magnetic field, 
may be essential to the formation of the observed structures. 

Mixing of cold gas may have a significant effect on the energy 
distribution in all cool core clusters. Of the nearby large cooling 
flows, it is only Abell 2029 which does not sh ow evidence for a 
large filamentary nebula (Johnstone et al ]|1987|) . In the Virgo clus- 
ter, where both Ha fi laments jSparks et al.l2004h and a weak shock 
jporman et alj|200^ are seen, the shock is outside the region of 
strong filamentary emission. Unlike Perseus, there is no evidence 
of deviation from the standard Rankine-Hugoniot relations for the 
shock in Virgo. 

Finally, we must consider the possibility that the pressure 
across the density jump is in fact continuous with an increase in the 
thermal gas pressure fraction in the "preshock" region. This could 
be the case if the medium outside the front has a larger component 
of ram, magnetic field or cosmic ray pressure, for example. 

The geometry of the density jump seems to disfavour an in- 
crease in ram pressure; because the shock is a continuous feature 
around the inner region it is difficult to maintain a bulk flow of the 
material in the whole outer region. 

The viability of the observed feature being produced by a 
change in the magnetic field strength across the density front de- 



S lends greatly on the field strength in the inner region. iTavlor et all 
2006h find a magnetic field of around 25 fiG based on rotation 
measure measurements of the very inner regions of the cluster core, 
although they suggest t his field is ass ociated largely with the fila- 
mentary nebula, whilst [Sanders et all (2005) estimate a magnetic 
field of around 0.1 /iG from non-thermal X-ray emission. In or- 
der to provide a continuous pressure across the shock, we require 
^outci- " dinner ^ 6x10'^ /iG^. If the inner field is around 25 fiG this 
is a factor of around 3 increase in the field strength. For a smaller 
outer field the required increase is correspondingly larger and for 
a 0.1 fiG iimer field, magnetic pressure balance across the density 
front would require an implausible increase in the field strength of 
a factor of ~ 1000. 

If there were a significant change in either the magnetic field 
or non-thermal particle population in the gas moving across the 
density front, we would expect to see a corresponding feature in the 
radio mini-halo emission. No-such feature can obviously be iden- 
tified in the VLA data of Sijbring ( 1993), although we cannot en- 
tirely rule out a model in which the non-thermal gas is the largest 
fraction by volume and the "shock" feature is associated with a 
change in the number density of clouds of thermal gas. Neverthe- 
less, the difficulties associated with these constant pressure models 
lead us to believe that an outward-propagating weak shock wave is 
the most likely explanation for the observed feature. 

In conclusion, the high pressure regions surrounding the inner 
radio bubbles in the Perseus cluster terminate in spherical shocks. 
The density jumps by 31 percent at the shock front, yet any imme- 
diate temperature jump is only 3 ± 6 percent. This is inconsistent 
with a simple adiabatic shock. We are unable to find a definitive so- 
lution to the problem. Several possibilities have been investigated, 
with turbulent post-shock gas rapidly mixing with cool material as- 
sociated with the optical filaments being promising. Such mixing, 
and the 3D geometry of the bubbles, can prevent the central gas 
from overheating, despite repeated bubbling over hundreds of mil- 
lions of years. The excess thermal energy contained in these high 
pressure region is about 3.5 times the total work done in expanding 
the bubbles adiabatically and the power provided by the shock is 
enough to significantly offset the cooling in Perseus. 
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